home *** CD-ROM | disk | FTP | other *** search
/ InfoMagic Internet Tools 1993 July / Internet Tools.iso / RockRidge / security / cops / cops_104 / src / crack-fcrypt.c < prev    next >
Encoding:
C/C++ Source or Header  |  1992-03-10  |  22.2 KB  |  833 lines
  1. /*
  2.  * This program is copyright Alec Muffett 1991 except for some portions of
  3.  * code in "crack-fcrypt.c" which are copyright Robert Baldwin, Icarus
  4.  * Sparry and Alec Muffett.  The author(s) disclaims all responsibility or
  5.  * liability with respect to it's usage or its effect upon hardware or
  6.  * computer systems, and maintain copyright as set out in the "LICENCE"
  7.  * document which accompanies distributions of Crack v4.0 and upwards. 
  8.  */
  9.  
  10. /*
  11.  * Misc defs for the fast password transform optimisations.
  12.  */
  13.  
  14. #include "crack.h"
  15.  
  16. /*
  17.  * Data segment gathered into one place, try to keep this stuff long aligned
  18.  * - AEM
  19.  */
  20.  
  21. static char iobuf[16];
  22. static obpb1 crypt_block[72];    /* 72 is next multiple of 8 bytes after 66 */
  23. static sbpb24 KS[32];
  24. static sbpb24 S0H[64], S1H[64], S2H[64], S3H[64];
  25. static sbpb24 S4H[64], S5H[64], S6H[64], S7H[64];
  26. static sbpb24 S0L[64], S1L[64], S2L[64], S3L[64];
  27. static sbpb24 S4L[64], S5L[64], S6L[64], S7L[64];
  28. static sbpb24 out96[4];
  29.  
  30. /*
  31.  * Start of the real thing
  32.  */
  33.  
  34. void
  35. fsetkey ()
  36. {
  37.     /*
  38.      * This used to be utterly horrendous. It still is, but it's much, much,
  39.      * smaller... (and quite a bit faster...) - AEM
  40.      */
  41.     static unsb KeyToKS[] =
  42.     {
  43.     9, 50, 33, 59, 48, 16, 32, 56, 1, 8, 18, 41, 2, 34, 25, 24,
  44.     43, 57, 58, 0, 35, 26, 17, 40, 21, 27, 38, 53, 36, 3, 46, 29,
  45.     4, 52, 22, 28, 60, 20, 37, 62, 14, 19, 44, 13, 12, 61, 54, 30,
  46.     1, 42, 25, 51, 40, 8, 24, 48, 58, 0, 10, 33, 59, 26, 17, 16,
  47.     35, 49, 50, 57, 56, 18, 9, 32, 13, 19, 30, 45, 28, 62, 38, 21,
  48.     27, 44, 14, 20, 52, 12, 29, 54, 6, 11, 36, 5, 4, 53, 46, 22,
  49.     50, 26, 9, 35, 24, 57, 8, 32, 42, 49, 59, 17, 43, 10, 1, 0,
  50.     48, 33, 34, 41, 40, 2, 58, 16, 60, 3, 14, 29, 12, 46, 22, 5,
  51.     11, 28, 61, 4, 36, 27, 13, 38, 53, 62, 20, 52, 19, 37, 30, 6,
  52.     34, 10, 58, 48, 8, 41, 57, 16, 26, 33, 43, 1, 56, 59, 50, 49,
  53.     32, 17, 18, 25, 24, 51, 42, 0, 44, 54, 61, 13, 27, 30, 6, 52,
  54.     62, 12, 45, 19, 20, 11, 60, 22, 37, 46, 4, 36, 3, 21, 14, 53,
  55.     18, 59, 42, 32, 57, 25, 41, 0, 10, 17, 56, 50, 40, 43, 34, 33,
  56.     16, 1, 2, 9, 8, 35, 26, 49, 28, 38, 45, 60, 11, 14, 53, 36,
  57.     46, 27, 29, 3, 4, 62, 44, 6, 21, 30, 19, 20, 54, 5, 61, 37,
  58.     2, 43, 26, 16, 41, 9, 25, 49, 59, 1, 40, 34, 24, 56, 18, 17,
  59.     0, 50, 51, 58, 57, 48, 10, 33, 12, 22, 29, 44, 62, 61, 37, 20,
  60.     30, 11, 13, 54, 19, 46, 28, 53, 5, 14, 3, 4, 38, 52, 45, 21,
  61.     51, 56, 10, 0, 25, 58, 9, 33, 43, 50, 24, 18, 8, 40, 2, 1,
  62.     49, 34, 35, 42, 41, 32, 59, 17, 27, 6, 13, 28, 46, 45, 21, 4,
  63.     14, 62, 60, 38, 3, 30, 12, 37, 52, 61, 54, 19, 22, 36, 29, 5,
  64.     35, 40, 59, 49, 9, 42, 58, 17, 56, 34, 8, 2, 57, 24, 51, 50,
  65.     33, 18, 48, 26, 25, 16, 43, 1, 11, 53, 60, 12, 30, 29, 5, 19,
  66.     61, 46, 44, 22, 54, 14, 27, 21, 36, 45, 38, 3, 6, 20, 13, 52,
  67.     56, 32, 51, 41, 1, 34, 50, 9, 48, 26, 0, 59, 49, 16, 43, 42,
  68.     25, 10, 40, 18, 17, 8, 35, 58, 3, 45, 52, 4, 22, 21, 60, 11,
  69.     53, 38, 36, 14, 46, 6, 19, 13, 28, 37, 30, 62, 61, 12, 5, 44,
  70.     40, 16, 35, 25, 50, 18, 34, 58, 32, 10, 49, 43, 33, 0, 56, 26,
  71.     9, 59, 24, 2, 1, 57, 48, 42, 54, 29, 36, 19, 6, 5, 44, 62,
  72.     37, 22, 20, 61, 30, 53, 3, 60, 12, 21, 14, 46, 45, 27, 52, 28,
  73.     24, 0, 48, 9, 34, 2, 18, 42, 16, 59, 33, 56, 17, 49, 40, 10,
  74.     58, 43, 8, 51, 50, 41, 32, 26, 38, 13, 20, 3, 53, 52, 28, 46,
  75.     21, 6, 4, 45, 14, 37, 54, 44, 27, 5, 61, 30, 29, 11, 36, 12,
  76.     8, 49, 32, 58, 18, 51, 2, 26, 0, 43, 17, 40, 1, 33, 24, 59,
  77.     42, 56, 57, 35, 34, 25, 16, 10, 22, 60, 4, 54, 37, 36, 12, 30,
  78.     5, 53, 19, 29, 61, 21, 38, 28, 11, 52, 45, 14, 13, 62, 20, 27,
  79.     57, 33, 16, 42, 2, 35, 51, 10, 49, 56, 1, 24, 50, 17, 8, 43,
  80.     26, 40, 41, 48, 18, 9, 0, 59, 6, 44, 19, 38, 21, 20, 27, 14,
  81.     52, 37, 3, 13, 45, 5, 22, 12, 62, 36, 29, 61, 60, 46, 4, 11,
  82.     41, 17, 0, 26, 51, 48, 35, 59, 33, 40, 50, 8, 34, 1, 57, 56,
  83.     10, 24, 25, 32, 2, 58, 49, 43, 53, 28, 3, 22, 5, 4, 11, 61,
  84.     36, 21, 54, 60, 29, 52, 6, 27, 46, 20, 13, 45, 44, 30, 19, 62,
  85.     25, 1, 49, 10, 35, 32, 48, 43, 17, 24, 34, 57, 18, 50, 41, 40,
  86.     59, 8, 9, 16, 51, 42, 33, 56, 37, 12, 54, 6, 52, 19, 62, 45,
  87.     20, 5, 38, 44, 13, 36, 53, 11, 30, 4, 60, 29, 28, 14, 3, 46,
  88.     17, 58, 41, 2, 56, 24, 40, 35, 9, 16, 26, 49, 10, 42, 33, 32,
  89.     51, 0, 1, 8, 43, 34, 25, 48, 29, 4, 46, 61, 44, 11, 54, 37,
  90.     12, 60, 30, 36, 5, 28, 45, 3, 22, 27, 52, 21, 20, 6, 62, 38
  91.     };
  92.  
  93.     reg int i;
  94.     reg int j;
  95.     reg unsigned long r;
  96.     reg unsb *k;
  97.  
  98.     k = KeyToKS;
  99.  
  100.     for (i = 0; i < 32; i++)
  101.     {
  102.     /* 16-bit tweaks suggested by cip_maz@fb6tcp.physik.uni-paderborn.de */
  103.     /* inlining speedup tweak suggested by tahorsley@csd.harris.com */
  104.     /* (strange addition compensates missing TF_TO_SIXBIT) */
  105.     r = (unsigned long) crypt_block[*(k++)];
  106.     r |= (unsigned long) crypt_block[*(k++)] << 1;
  107.     r |= (unsigned long) crypt_block[*(k++)] << 2;
  108.     r |= (unsigned long) crypt_block[*(k++)] << 3;
  109.     r |= (unsigned long) crypt_block[*(k++)] << 4;
  110.     r |= (unsigned long) crypt_block[*(k++)] << 5;
  111.     r |= (unsigned long) crypt_block[*(k++)] << (2 + 6);
  112.     r |= (unsigned long) crypt_block[*(k++)] << (2 + 7);
  113.     r |= (unsigned long) crypt_block[*(k++)] << (2 + 8);
  114.     r |= (unsigned long) crypt_block[*(k++)] << (2 + 9);
  115.     r |= (unsigned long) crypt_block[*(k++)] << (2 + 10);
  116.     r |= (unsigned long) crypt_block[*(k++)] << (2 + 11);
  117.     r |= (unsigned long) crypt_block[*(k++)] << (4 + 12);
  118.     r |= (unsigned long) crypt_block[*(k++)] << (4 + 13);
  119.     r |= (unsigned long) crypt_block[*(k++)] << (4 + 14);
  120.     r |= (unsigned long) crypt_block[*(k++)] << (4 + 15);
  121.     r |= (unsigned long) crypt_block[*(k++)] << (4 + 16);
  122.     r |= (unsigned long) crypt_block[*(k++)] << (4 + 17);
  123.     r |= (unsigned long) crypt_block[*(k++)] << (6 + 18);
  124.     r |= (unsigned long) crypt_block[*(k++)] << (6 + 19);
  125.     r |= (unsigned long) crypt_block[*(k++)] << (6 + 20);
  126.     r |= (unsigned long) crypt_block[*(k++)] << (6 + 21);
  127.     r |= (unsigned long) crypt_block[*(k++)] << (6 + 22);
  128.     r |= (unsigned long) crypt_block[*(k++)] << (6 + 23);
  129.     KS[i] = r;
  130.     }
  131. }
  132.  
  133. void
  134. XForm (saltvalue)
  135.     sbpb24 saltvalue;
  136. {
  137. #ifdef BIG_ENDIAN        /* Icarus Sparry, Bath - mod AEM */
  138. #define STEP --
  139. #define START &sdata.c[7]
  140. #define Dl sdata.b[1]
  141. #define Dh sdata.b[0]
  142. #else                /* LITTLE_ENDIAN */
  143. #define STEP ++
  144. #define START &sdata.c[0]
  145. #define Dl sdata.b[0]
  146. #define Dh sdata.b[1]
  147. #endif
  148.     union SDATA sdata;
  149.     reg sbpb24 Rl;
  150.     reg sbpb24 Rh;
  151.     reg sbpb24 Ll;
  152.     reg sbpb24 Lh;
  153.     reg sbpb6 *dp;
  154.     int loop;
  155.     int kloop;
  156.     sbpb24 *kp;
  157.     reg sbpb24 k;
  158. #ifdef FDES_8BYTE
  159.     reg sbpb24 tmpi;
  160. #endif    /* FDES_8BYTE */
  161.  
  162.     Ll = Lh = Rl = Rh = 0;
  163.  
  164.     for (loop = 25; loop--; /* nothing */ )
  165.     {
  166.     kp = KS;
  167.     for (kloop = 8; kloop--; /* nothing */ )
  168.     {
  169.         k = (Rl ^ Rh) & saltvalue;
  170. #ifndef FDES_8BYTE
  171.         Dl = (k ^ Rl ^ *kp++) << SIZEFIX;
  172.         Dh = (k ^ Rh ^ *kp++) << SIZEFIX;
  173. #else
  174.         /* hack to make things work better - matthew kaufman */
  175.         /* I haven't tried any of this - I don't have a cray... AEM */
  176.         tmpi = (k ^ Rl ^ *kp++);
  177.         sdata.c[3] = (tmpi >> 24) & 0x00ff;
  178.         sdata.c[2] = (tmpi >> 16) & 0x00ff;
  179.         sdata.c[1] = (tmpi >> 8) & 0x00ff;
  180.         sdata.c[0] = (tmpi) & 0x00ff;
  181.         tmpi = (k ^ Rh ^ *kp++);
  182.         sdata.c[7] = (tmpi >> 24) & 0x00ff;
  183.         sdata.c[6] = (tmpi >> 16) & 0x00ff;
  184.         sdata.c[5] = (tmpi >> 8) & 0x00ff;
  185.         sdata.c[4] = (tmpi) & 0x00ff;
  186. #endif    /* FDES_8BYTE */
  187.  
  188.         dp = START;
  189.         Lh ^= INDIRECT (S0H, *dp);
  190.         Ll ^= INDIRECT (S0L, *dp STEP);
  191.         Lh ^= INDIRECT (S1H, *dp);
  192.         Ll ^= INDIRECT (S1L, *dp STEP);
  193.         Lh ^= INDIRECT (S2H, *dp);
  194.         Ll ^= INDIRECT (S2L, *dp STEP);
  195.         Lh ^= INDIRECT (S3H, *dp);
  196.         Ll ^= INDIRECT (S3L, *dp STEP);
  197.         Lh ^= INDIRECT (S4H, *dp);
  198.         Ll ^= INDIRECT (S4L, *dp STEP);
  199.         Lh ^= INDIRECT (S5H, *dp);
  200.         Ll ^= INDIRECT (S5L, *dp STEP);
  201.         Lh ^= INDIRECT (S6H, *dp);
  202.         Ll ^= INDIRECT (S6L, *dp STEP);
  203.         Lh ^= INDIRECT (S7H, *dp);
  204.         Ll ^= INDIRECT (S7L, *dp STEP);
  205.  
  206.         k = (Ll ^ Lh) & saltvalue;
  207. #ifndef FDES_8BYTE
  208.         Dl = (k ^ Ll ^ *kp++) << SIZEFIX;
  209.         Dh = (k ^ Lh ^ *kp++) << SIZEFIX;
  210. #else
  211.         tmpi = (k ^ Ll ^ *kp++);
  212.         sdata.c[3] = (tmpi >> 24) & 0x00ff;
  213.         sdata.c[2] = (tmpi >> 16) & 0x00ff;
  214.         sdata.c[1] = (tmpi >> 8) & 0x00ff;
  215.         sdata.c[0] = (tmpi) & 0x00ff;
  216.         tmpi = (k ^ Lh ^ *kp++);
  217.         sdata.c[7] = (tmpi >> 24) & 0x00ff;
  218.         sdata.c[6] = (tmpi >> 16) & 0x00ff;
  219.         sdata.c[5] = (tmpi >> 8) & 0x00ff;
  220.         sdata.c[4] = (tmpi) & 0x00ff;
  221. #endif    /* FDES_8BYTE */
  222.  
  223.         dp = START;
  224.         Rh ^= INDIRECT (S0H, *dp);
  225.         Rl ^= INDIRECT (S0L, *dp STEP);
  226.         Rh ^= INDIRECT (S1H, *dp);
  227.         Rl ^= INDIRECT (S1L, *dp STEP);
  228.         Rh ^= INDIRECT (S2H, *dp);
  229.         Rl ^= INDIRECT (S2L, *dp STEP);
  230.         Rh ^= INDIRECT (S3H, *dp);
  231.         Rl ^= INDIRECT (S3L, *dp STEP);
  232.         Rh ^= INDIRECT (S4H, *dp);
  233.         Rl ^= INDIRECT (S4L, *dp STEP);
  234.         Rh ^= INDIRECT (S5H, *dp);
  235.         Rl ^= INDIRECT (S5L, *dp STEP);
  236.         Rh ^= INDIRECT (S6H, *dp);
  237.         Rl ^= INDIRECT (S6L, *dp STEP);
  238.         Rh ^= INDIRECT (S7H, *dp);
  239.         Rl ^= INDIRECT (S7L, *dp STEP);
  240.     }
  241.  
  242.     Ll ^= Rl;
  243.     Lh ^= Rh;
  244.     Rl ^= Ll;
  245.     Rh ^= Lh;
  246.     Ll ^= Rl;
  247.     Lh ^= Rh;
  248.     }
  249.  
  250.     /*
  251.      * for reasons that I cannot explain, if I insert the contents of the
  252.      * UnXForm function right HERE, making the tweaks as necessary to avoid
  253.      * using out96[] to pass data, I LOSE 30% of my speed.  I don't know why.
  254.      * Hence, I continue to use out96[]...
  255.      */
  256.     {
  257.     reg sbpb24 *qp;
  258.     qp = out96;
  259.     *qp++ = Ll;
  260.     *qp++ = Lh;
  261.     *qp++ = Rl;
  262.     *qp++ = Rh;
  263.     }
  264. }
  265.  
  266. void
  267. UnXForm ()
  268. {
  269.     reg sbpb24 Rl;
  270.     reg sbpb24 Rh;
  271.     reg sbpb24 Ll;
  272.     reg sbpb24 Lh;
  273.     reg obpb1 *ptr;
  274.     reg long int mask;
  275.     register long int *lip;
  276.  
  277.     Ll = SIXBIT_TO_TF (out96[0]);
  278.     Lh = SIXBIT_TO_TF (out96[1]);
  279.     Rl = SIXBIT_TO_TF (out96[2]);
  280.     Rh = SIXBIT_TO_TF (out96[3]);
  281.  
  282. #ifdef BUILTIN_CLEAR
  283.     lip = (long int *) crypt_block;
  284.     for (mask = (sizeof (crypt_block) / sizeof (long int)); mask--; /* - */ )
  285.     {
  286.     *(lip++) = 0L;
  287.     }
  288. #else                /* BUILTIN_CLEAR */
  289.     bzero (crypt_block, 66);
  290. #endif                /* BUILTIN_CLEAR */
  291.  
  292.     ptr = crypt_block;
  293.     mask = 0x000400L;
  294.     if (Rl & mask)
  295.     *ptr = 0x01;
  296.     ptr++;
  297.     if (Ll & mask)
  298.     *ptr = 0x01;
  299.     ptr++;
  300.     mask = 0x400000L;
  301.     if (Rl & mask)
  302.     *ptr = 0x01;
  303.     ptr++;
  304.     if (Ll & mask)
  305.     *ptr = 0x01;
  306.     ptr++;
  307.     mask = 0x000400L;
  308.     if (Rh & mask)
  309.     *ptr = 0x01;
  310.     ptr++;
  311.     if (Lh & mask)
  312.     *ptr = 0x01;
  313.     ptr++;
  314.     mask = 0x400000L;
  315.     if (Rh & mask)
  316.     *ptr = 0x01;
  317.     ptr++;
  318.     if (Lh & mask)
  319.     *ptr = 0x01;
  320.     ptr++;
  321.     mask = 0x000200L;
  322.     if (Rl & mask)
  323.     *ptr = 0x01;
  324.     ptr++;
  325.     if (Ll & mask)
  326.     *ptr = 0x01;
  327.     ptr++;
  328.     mask = 0x200000L;
  329.     if (Rl & mask)
  330.     *ptr = 0x01;
  331.     ptr++;
  332.     if (Ll & mask)
  333.     *ptr = 0x01;
  334.     ptr++;
  335.     mask = 0x000200L;
  336.     if (Rh & mask)
  337.     *ptr = 0x01;
  338.     ptr++;
  339.     if (Lh & mask)
  340.     *ptr = 0x01;
  341.     ptr++;
  342.     mask = 0x200000L;
  343.     if (Rh & mask)
  344.     *ptr = 0x01;
  345.     ptr++;
  346.     if (Lh & mask)
  347.     *ptr = 0x01;
  348.     ptr++;
  349.     mask = 0x000100L;
  350.     if (Rl & mask)
  351.     *ptr = 0x01;
  352.     ptr++;
  353.     if (Ll & mask)
  354.     *ptr = 0x01;
  355.     ptr++;
  356.     mask = 0x100000L;
  357.     if (Rl & mask)
  358.     *ptr = 0x01;
  359.     ptr++;
  360.     if (Ll & mask)
  361.     *ptr = 0x01;
  362.     ptr++;
  363.     mask = 0x000100L;
  364.     if (Rh & mask)
  365.     *ptr = 0x01;
  366.     ptr++;
  367.     if (Lh & mask)
  368.     *ptr = 0x01;
  369.     ptr++;
  370.     mask = 0x100000L;
  371.     if (Rh & mask)
  372.     *ptr = 0x01;
  373.     ptr++;
  374.     if (Lh & mask)
  375.     *ptr = 0x01;
  376.     ptr++;
  377.     mask = 0x000080L;
  378.     if (Rl & mask)
  379.     *ptr = 0x01;
  380.     ptr++;
  381.     if (Ll & mask)
  382.     *ptr = 0x01;
  383.     ptr++;
  384.     mask = 0x080000L;
  385.     if (Rl & mask)
  386.     *ptr = 0x01;
  387.     ptr++;
  388.     if (Ll & mask)
  389.     *ptr = 0x01;
  390.     ptr++;
  391.     mask = 0x000080L;
  392.     if (Rh & mask)
  393.     *ptr = 0x01;
  394.     ptr++;
  395.     if (Lh & mask)
  396.     *ptr = 0x01;
  397.     ptr++;
  398.     mask = 0x080000L;
  399.     if (Rh & mask)
  400.     *ptr = 0x01;
  401.     ptr++;
  402.     if (Lh & mask)
  403.     *ptr = 0x01;
  404.     ptr++;
  405.     mask = 0x000010L;
  406.     if (Rl & mask)
  407.     *ptr = 0x01;
  408.     ptr++;
  409.     if (Ll & mask)
  410.     *ptr = 0x01;
  411.     ptr++;
  412.     mask = 0x010000L;
  413.     if (Rl & mask)
  414.     *ptr = 0x01;
  415.     ptr++;
  416.     if (Ll & mask)
  417.     *ptr = 0x01;
  418.     ptr++;
  419.     mask = 0x000010L;
  420.     if (Rh & mask)
  421.     *ptr = 0x01;
  422.     ptr++;
  423.     if (Lh & mask)
  424.     *ptr = 0x01;
  425.     ptr++;
  426.     mask = 0x010000L;
  427.     if (Rh & mask)
  428.     *ptr = 0x01;
  429.     ptr++;
  430.     if (Lh & mask)
  431.     *ptr = 0x01;
  432.     ptr++;
  433.     mask = 0x000008L;
  434.     if (Rl & mask)
  435.     *ptr = 0x01;
  436.     ptr++;
  437.     if (Ll & mask)
  438.     *ptr = 0x01;
  439.     ptr++;
  440.     mask = 0x008000L;
  441.     if (Rl & mask)
  442.     *ptr = 0x01;
  443.     ptr++;
  444.     if (Ll & mask)
  445.     *ptr = 0x01;
  446.     ptr++;
  447.     mask = 0x000008L;
  448.     if (Rh & mask)
  449.     *ptr = 0x01;
  450.     ptr++;
  451.     if (Lh & mask)
  452.     *ptr = 0x01;
  453.     ptr++;
  454.     mask = 0x008000L;
  455.     if (Rh & mask)
  456.     *ptr = 0x01;
  457.     ptr++;
  458.     if (Lh & mask)
  459.     *ptr = 0x01;
  460.     ptr++;
  461.     mask = 0x000004L;
  462.     if (Rl & mask)
  463.     *ptr = 0x01;
  464.     ptr++;
  465.     if (Ll & mask)
  466.     *ptr = 0x01;
  467.     ptr++;
  468.     mask = 0x004000L;
  469.     if (Rl & mask)
  470.     *ptr = 0x01;
  471.     ptr++;
  472.     if (Ll & mask)
  473.     *ptr = 0x01;
  474.     ptr++;
  475.     mask = 0x000004L;
  476.     if (Rh & mask)
  477.     *ptr = 0x01;
  478.     ptr++;
  479.     if (Lh & mask)
  480.     *ptr = 0x01;
  481.     ptr++;
  482.     mask = 0x004000L;
  483.     if (Rh & mask)
  484.     *ptr = 0x01;
  485.     ptr++;
  486.     if (Lh & mask)
  487.     *ptr = 0x01;
  488.     ptr++;
  489.     mask = 0x000002L;
  490.     if (Rl & mask)
  491.     *ptr = 0x01;
  492.     ptr++;
  493.     if (Ll & mask)
  494.     *ptr = 0x01;
  495.     ptr++;
  496.     mask = 0x002000L;
  497.     if (Rl & mask)
  498.     *ptr = 0x01;
  499.     ptr++;
  500.     if (Ll & mask)
  501.     *ptr = 0x01;
  502.     ptr++;
  503.     mask = 0x000002L;
  504.     if (Rh & mask)
  505.     *ptr = 0x01;
  506.     ptr++;
  507.     if (Lh & mask)
  508.     *ptr = 0x01;
  509.     ptr++;
  510.     mask = 0x002000L;
  511.     if (Rh & mask)
  512.     *ptr = 0x01;
  513.     ptr++;
  514.     if (Lh & mask)
  515.     *ptr = 0x01;
  516.     ptr++;
  517. }
  518.  
  519. char *
  520. fcrypt (pw, salt)
  521.     char *pw;
  522.     char *salt;
  523. {
  524.     /* Table lookups for salts reduce fcrypt() overhead dramatically */
  525.     static sbpb24 salt0[] =
  526.     {
  527.     18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31,
  528.     32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47,
  529.     48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63,
  530.     0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 5, 6, 7, 8, 9, 10, 11, 12,
  531.     13, 14, 15, 16, 17,
  532.  
  533.     18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31,
  534.     32, 33, 34, 35, 36, 37, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41,
  535.     42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57,
  536.     58, 59, 60, 61, 62, 63, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,
  537.     12, 13, 14, 15, 16, 17,
  538.  
  539.     18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31,
  540.     32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47,
  541.     48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63,
  542.     0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17,
  543.  
  544.     18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31,
  545.     32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47,
  546.     48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63,
  547.     0, 1, 2, 3, 4
  548.     };
  549.     static sbpb24 salt1[] =
  550.     {
  551.     1152, 1216, 1280, 1344, 1408, 1472, 1536, 1600, 1664,
  552.     1728, 1792, 1856, 1920, 1984, 2048, 2112, 2176, 2240, 2304,
  553.     2368, 2432, 2496, 2560, 2624, 2688, 2752, 2816, 2880, 2944,
  554.     3008, 3072, 3136, 3200, 3264, 3328, 3392, 3456, 3520, 3584,
  555.     3648, 3712, 3776, 3840, 3904, 3968, 4032, 0, 64, 128, 192, 256,
  556.     320, 384, 448, 512, 576, 640, 704, 320, 384, 448, 512, 576, 640,
  557.     704, 768, 832, 896, 960, 1024, 1088,
  558.  
  559.     1152, 1216, 1280, 1344, 1408, 1472, 1536, 1600, 1664,
  560.     1728, 1792, 1856, 1920, 1984, 2048, 2112, 2176, 2240, 2304,
  561.     2368, 2048, 2112, 2176, 2240, 2304, 2368, 2432, 2496, 2560,
  562.     2624, 2688, 2752, 2816, 2880, 2944, 3008, 3072, 3136, 3200,
  563.     3264, 3328, 3392, 3456, 3520, 3584, 3648, 3712, 3776, 3840,
  564.     3904, 3968, 4032, 0, 64, 128, 192, 256, 320, 384, 448, 512, 576,
  565.     640, 704, 768, 832, 896, 960, 1024, 1088,
  566.  
  567.     1152, 1216, 1280, 1344, 1408, 1472, 1536, 1600, 1664,
  568.     1728, 1792, 1856, 1920, 1984, 2048, 2112, 2176, 2240, 2304,
  569.     2368, 2432, 2496, 2560, 2624, 2688, 2752, 2816, 2880, 2944,
  570.     3008, 3072, 3136, 3200, 3264, 3328, 3392, 3456, 3520, 3584,
  571.     3648, 3712, 3776, 3840, 3904, 3968, 4032, 0, 64, 128, 192, 256,
  572.     320, 384, 448, 512, 576, 640, 704, 768, 832, 896, 960, 1024,
  573.     1088,
  574.  
  575.     1152, 1216, 1280, 1344, 1408, 1472, 1536, 1600, 1664,
  576.     1728, 1792, 1856, 1920, 1984, 2048, 2112, 2176, 2240, 2304,
  577.     2368, 2432, 2496, 2560, 2624, 2688, 2752, 2816, 2880, 2944,
  578.     3008, 3072, 3136, 3200, 3264, 3328, 3392, 3456, 3520, 3584,
  579.     3648, 3712, 3776, 3840, 3904, 3968, 4032, 0, 64, 128, 192, 256
  580.     };
  581.  
  582.     /* final perutation desalting */
  583.     static obpb1 final[] =
  584.     {
  585.     46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 65, 66,
  586.     67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82,
  587.     83, 84, 85, 86, 87, 88, 89, 90, 97, 98, 99, 100, 101, 102, 103,
  588.     104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116,
  589.     117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129,
  590.     130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142,
  591.     143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155,
  592.     156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168,
  593.     169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181,
  594.     182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194,
  595.     195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 207,
  596.     208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220,
  597.     221, 222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233,
  598.     234, 235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245, 246,
  599.     247, 248, 249, 250, 251, 252, 253, 254, 255,
  600.         /* Truncate overflow bits at 256 */
  601.     0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15,
  602.     16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31,
  603.     32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47,
  604.     48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58
  605.     };
  606.  
  607.     reg int i, j, k;
  608.     reg long int *lip;
  609.     sbpb24 saltvalue;
  610.  
  611. #ifdef BUILTIN_CLEAR
  612.     lip = (long int *) crypt_block;
  613.     for (i = (sizeof (crypt_block) / sizeof (long int)); i--; /* - */ )
  614.     {
  615.     *(lip++) = 0L;
  616.     }
  617. #else                /* BUILTIN_CLEAR */
  618.     bzero (crypt_block, 66);
  619. #endif                /* BUILTIN_CLEAR */
  620.  
  621.     for (i = 0; (k = *pw) && i < 64; pw++)
  622.     {
  623.     crypt_block[i++] = (k >> 6) & 01;
  624.     crypt_block[i++] = (k >> 5) & 01;
  625.     crypt_block[i++] = (k >> 4) & 01;
  626.     crypt_block[i++] = (k >> 3) & 01;
  627.     crypt_block[i++] = (k >> 2) & 01;
  628.     crypt_block[i++] = (k >> 1) & 01;
  629.     crypt_block[i++] = (k >> 0) & 01;
  630.     i++;            /* have to skip one here (parity bit) */
  631.     }
  632.  
  633.     fsetkey ( /* crypt_block */ );
  634.  
  635. #ifdef BUILTIN_CLEAR
  636.     lip = (long int *) crypt_block;
  637.     for (i = (sizeof (crypt_block) / sizeof (long int)); i--; /* - */ )
  638.     {
  639.     *(lip++) = 0L;
  640.     }
  641. #else                /* BUILTIN_CLEAR */
  642.     bzero (crypt_block, 66);
  643. #endif                /* BUILTIN_CLEAR */
  644.  
  645.     iobuf[0] = salt[0];
  646.     iobuf[1] = salt[1];
  647.  
  648.     saltvalue = salt0[iobuf[0]] | salt1[iobuf[1]];
  649.     saltvalue = TF_TO_SIXBIT (saltvalue);
  650.  
  651.     XForm (saltvalue);
  652.     UnXForm ();
  653.  
  654.     for (i = 0; i < 11; i++)
  655.     {
  656.     k = 0;
  657.  
  658.     for (j = 0; j < 6; j++)
  659.     {
  660.         k = (k << 1) | crypt_block[SIX_TIMES (i) + j];
  661.     }
  662.     iobuf[i + 2] = final[k];
  663.     }
  664.  
  665.     iobuf[i + 2] = 0;
  666.  
  667.     if (iobuf[1] == 0)
  668.     {
  669.     iobuf[1] = iobuf[0];
  670.     }
  671.     return (iobuf);
  672. }
  673. /********* INITIALISATION ROUTINES *********/
  674.  
  675. fbpb4
  676. lookupS (tableno, t6bits)
  677.     unsl tableno;
  678.     sbpb6R t6bits;
  679. {
  680.     static fbpb4R S[8][64] =
  681.     {
  682.     14, 4, 13, 1, 2, 15, 11, 8, 3, 10, 6, 12, 5, 9, 0, 7,
  683.     0, 15, 7, 4, 14, 2, 13, 1, 10, 6, 12, 11, 9, 5, 3, 8,
  684.     4, 1, 14, 8, 13, 6, 2, 11, 15, 12, 9, 7, 3, 10, 5, 0,
  685.     15, 12, 8, 2, 4, 9, 1, 7, 5, 11, 3, 14, 10, 0, 6, 13,
  686.  
  687.     15, 1, 8, 14, 6, 11, 3, 4, 9, 7, 2, 13, 12, 0, 5, 10,
  688.     3, 13, 4, 7, 15, 2, 8, 14, 12, 0, 1, 10, 6, 9, 11, 5,
  689.     0, 14, 7, 11, 10, 4, 13, 1, 5, 8, 12, 6, 9, 3, 2, 15,
  690.     13, 8, 10, 1, 3, 15, 4, 2, 11, 6, 7, 12, 0, 5, 14, 9,
  691.  
  692.     10, 0, 9, 14, 6, 3, 15, 5, 1, 13, 12, 7, 11, 4, 2, 8,
  693.     13, 7, 0, 9, 3, 4, 6, 10, 2, 8, 5, 14, 12, 11, 15, 1,
  694.     13, 6, 4, 9, 8, 15, 3, 0, 11, 1, 2, 12, 5, 10, 14, 7,
  695.     1, 10, 13, 0, 6, 9, 8, 7, 4, 15, 14, 3, 11, 5, 2, 12,
  696.  
  697.     7, 13, 14, 3, 0, 6, 9, 10, 1, 2, 8, 5, 11, 12, 4, 15,
  698.     13, 8, 11, 5, 6, 15, 0, 3, 4, 7, 2, 12, 1, 10, 14, 9,
  699.     10, 6, 9, 0, 12, 11, 7, 13, 15, 1, 3, 14, 5, 2, 8, 4,
  700.     3, 15, 0, 6, 10, 1, 13, 8, 9, 4, 5, 11, 12, 7, 2, 14,
  701.  
  702.     2, 12, 4, 1, 7, 10, 11, 6, 8, 5, 3, 15, 13, 0, 14, 9,
  703.     14, 11, 2, 12, 4, 7, 13, 1, 5, 0, 15, 10, 3, 9, 8, 6,
  704.     4, 2, 1, 11, 10, 13, 7, 8, 15, 9, 12, 5, 6, 3, 0, 14,
  705.     11, 8, 12, 7, 1, 14, 2, 13, 6, 15, 0, 9, 10, 4, 5, 3,
  706.  
  707.     12, 1, 10, 15, 9, 2, 6, 8, 0, 13, 3, 4, 14, 7, 5, 11,
  708.     10, 15, 4, 2, 7, 12, 9, 5, 6, 1, 13, 14, 0, 11, 3, 8,
  709.     9, 14, 15, 5, 2, 8, 12, 3, 7, 0, 4, 10, 1, 13, 11, 6,
  710.     4, 3, 2, 12, 9, 5, 15, 10, 11, 14, 1, 7, 6, 0, 8, 13,
  711.  
  712.     4, 11, 2, 14, 15, 0, 8, 13, 3, 12, 9, 7, 5, 10, 6, 1,
  713.     13, 0, 11, 7, 4, 9, 1, 10, 14, 3, 5, 12, 2, 15, 8, 6,
  714.     1, 4, 11, 13, 12, 3, 7, 14, 10, 15, 6, 8, 0, 5, 9, 2,
  715.     6, 11, 13, 8, 1, 4, 10, 7, 9, 5, 0, 15, 14, 2, 3, 12,
  716.  
  717.     13, 2, 8, 4, 6, 15, 11, 1, 10, 9, 3, 14, 5, 0, 12, 7,
  718.     1, 15, 13, 8, 10, 3, 7, 4, 12, 5, 6, 11, 0, 14, 9, 2,
  719.     7, 11, 4, 1, 9, 12, 14, 2, 0, 6, 10, 13, 15, 3, 5, 8,
  720.     2, 1, 14, 7, 4, 10, 8, 13, 15, 12, 9, 0, 3, 5, 6, 11,
  721.     };
  722.     sbpb6 fixed6bits;
  723.     fbpb4R r;
  724.     fbpb4 fixedr;
  725.  
  726.     fixed6bits = (((t6bits >> 0) & 01) << 5) +
  727.     (((t6bits >> 1) & 01) << 3) +
  728.     (((t6bits >> 2) & 01) << 2) +
  729.     (((t6bits >> 3) & 01) << 1) +
  730.     (((t6bits >> 4) & 01) << 0) +
  731.     (((t6bits >> 5) & 01) << 4);
  732.  
  733.     r = S[tableno][fixed6bits];
  734.  
  735.     fixedr = (((r >> 3) & 01) << 0) +
  736.     (((r >> 2) & 01) << 1) +
  737.     (((r >> 1) & 01) << 2) +
  738.     (((r >> 0) & 01) << 3);
  739.  
  740.     return (fixedr);
  741. }
  742.  
  743. void
  744. init (tableno, lowptr, highptr)
  745.     unsl tableno;
  746.     sbpb24 *lowptr, *highptr;
  747. {
  748.  
  749.     static unsb P[] =
  750.     {
  751.     15, 6, 19, 20,
  752.     28, 11, 27, 16,
  753.     0, 14, 22, 25,
  754.     4, 17, 30, 9,
  755.     1, 7, 23, 13,
  756.     31, 26, 2, 8,
  757.     18, 12, 29, 5,
  758.     21, 10, 3, 24,
  759.     };
  760.  
  761.     static unsb E[] =
  762.     {
  763.     31, 0, 1, 2, 3, 4,
  764.     3, 4, 5, 6, 7, 8,
  765.     7, 8, 9, 10, 11, 12,
  766.     11, 12, 13, 14, 15, 16,
  767.     15, 16, 17, 18, 19, 20,
  768.     19, 20, 21, 22, 23, 24,
  769.     23, 24, 25, 26, 27, 28,
  770.     27, 28, 29, 30, 31, 0,
  771.     };
  772.  
  773.     static obpb1 tmp32[32];
  774.     static obpb1 tmpP32[32];
  775.     static obpb1 tmpE[48];
  776.  
  777.     int j, k, i;
  778.     int tablenoX4;
  779.     reg sbpb24 spare24;
  780.  
  781.     tablenoX4 = tableno * 4;
  782.  
  783.     for (j = 0; j < 64; j++)
  784.     {
  785.     k = lookupS (tableno, j);
  786.  
  787.     for (i = 0; i < 32; i++)
  788.     {
  789.         tmp32[i] = 0;
  790.     }
  791.     for (i = 0; i < 4; i++)
  792.     {
  793.         tmp32[tablenoX4 + i] = (k >> i) & 01;
  794.     }
  795.     for (i = 0; i < 32; i++)
  796.     {
  797.         tmpP32[i] = tmp32[P[i]];
  798.     }
  799.     for (i = 0; i < 48; i++)
  800.     {
  801.         tmpE[i] = tmpP32[E[i]];
  802.     }
  803.  
  804.     lowptr[j] = 0;
  805.     highptr[j] = 0;
  806.  
  807.     for (i = 0; i < 24; i++)
  808.     {
  809.         lowptr[j] |= (unsigned long) tmpE[i] << i;
  810.     }
  811.     for (k = 0, i = 24; i < 48; i++, k++)
  812.     {
  813.         highptr[j] |= (unsigned long) tmpE[i] << k;
  814.     }
  815.  
  816.     spare24 = lowptr[j];    /* to allow for macro expansion */
  817.     lowptr[j] = TF_TO_SIXBIT (spare24);
  818.     spare24 = highptr[j];    /* to allow for macro expansion */
  819.     highptr[j] = TF_TO_SIXBIT (spare24);
  820.     }
  821. }
  822. init_des ()
  823. {
  824.     init (0L, S0L, S0H);
  825.     init (1L, S1L, S1H);
  826.     init (2L, S2L, S2H);
  827.     init (3L, S3L, S3H);
  828.     init (4L, S4L, S4H);
  829.     init (5L, S5L, S5H);
  830.     init (6L, S6L, S6H);
  831.     init (7L, S7L, S7H);
  832. }
  833.